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vi.  Overview 

This  report  provides  a  preliminary  summary  of  a 
geochemical  data  contribution  to  methane  hydrate 
research  and  exploration  on  the  Hikurangi  Margin,  off 
the  northeastern  coast  of  New  Zealand  from  June  20,  to 
July  3,  2006.  Geochemical  porewater  profiles  taken 
from  shallow  piston  cores  and  vertical  fluid  migration 
measured  with  heatflow  probing  were  compared  with 
seismic  summaries  of  potential  deep  sediment  hydrates 
deposits.  Recent  geochemical  and  geophysical 
investigation  of  hydrates  on  the  Cascadia  Margin,  mid- 
Chilean  coast,  and  Texas-Louisiana  Shelf  demonstrated 
a  need  to  combine  these  data  for  prediction  of  deep 
sediment  hydrate  deposits.  In  addition  to  methane 
hydrate  exploration  basic  sediment  and  water  column 
research  topics  were  addressed  during  this  cruise. 
General  research  goals  for  this  expedition  include: 

a) .  Refine  geophysical,  geochemical  and 
microbiological  technologies  for  prospecting  hydrate 
distribution  and  content. 

b) .  Contribute  to  establishing  high-priority 
geographical  regions  of  prospective  interest,  in  terms  of 
methane  volume  estimates. 

c) .  Prediction  of  environmental  effects  and  geologic 
risks  at  the  continental  margin  associated  to  the  natural 
resource  occurrence  and  resource  exploitation. 

d) .  Expand  understanding  of  the  biogeochemical 
parameters  and  associated  microbial  community 
diversity  in  shallow  sediments  that  influence  the 
porewater  sulfate  gradient  observed  through  anaerobic 
oxidation  of  methane. 

This  expedition  combines  several  ocean  research 
topics.  On  a  global  perspective,  methane  hydrates  are 
being  survey  by  many  nations  as  a  new  energy  source. 
Exploration  through  different  regions  of  the  world 
ocean  contributes  to  the  estimate  of  this  energy  resource 
relative  to  the  total  fossil  fuel  reserve.  Also  with 
recoginition  of  the  potential  coastal  sediment  hydrate 
loading  through  the  world  ocean  a  more  thorough 
understanding  of  ocean  carbon  and  climate  cycles 
requires  consideration  of  methane. 

This  project  continues  NRL  and  ONR-Global  goals 
for  international  collaboration  on  methane  hydrate 


research.  The  data  support  development  of  the 
NRL  program  on  methane  hydrate  influence  on 
sediment  geoacoustic  and  geotechnical  properties 
in  various  coastal  regions.  In  addition  the  program 
extends  the  current  data  base  from  recent  work  in 
coastal  regions  around  western  Canada  (Cascadia 
Margin),  the  Texas-Louisiana  Shelf,  Blake  Ridge, 
and  mid  Chilean  coast  which  show  diverse 
seismic  profiles  that  result  from  differences  in  the 
sediment  morphology  and  hydrate  gas  sources. 
Finally,  this  international  collaboration  shares  the 
high  cost  of  ship  time  and  development  of  new 
field  technologies. 

I.  Introduction 

Continental  margins  throughout  the  world  are 
rich  in  sediment  methane  hydrates.  The  amount 
of  methane  gas  trapped  in  hydrates  is  estimated  to 
be  significant  relative  to  the  energy  of  the  total 
conventional  fossil  fuel  reserves.  The  NRL 
methane  hydrate  research  focus  is  to  develop 
improved  models  for  study  of  the  creation  and 
dissociation  of  natural  gas  hydrates  and  to 
quantify  the  impact  of  these  processes  on  the 
geophysical  and  geotechnical  properties  of  marine 
sediments.  In  addition,  a  large  global  volume  of 
methane  promotes  interest  in  a  new  energy 
source,  coastal  stability,  alteration  of  global 
economy  and  understandings  of  time  scale  and 
steps  in  climate  change. 

Alteration  of  hydrate  stability  in  coastal 
sediments  can  impact  drilling  platform  stability 
and  contribute  to  environmental  hazards. 
Submarine  geo-hazards  (such  as  sediment 
instabilities  and  slope  failures  on  local  and 
regional  scales,  leading  to  debris  flows,  slumps, 
slides,  and  possible  tsunamis)  caused  by  gas- 
hydrate  dissociation  are  of  immediate  and 
increasing  importance  as  humankind  moves  to 
exploit  seabed  resources  in  ever-deepening  waters 
of  coastal  oceans. 

On  a  natural  global  cycle  the  vulnerability  of 
gas  hydrate  to  temperature  and  sea  level  change 
varies  the  stability  of  deep-water  oceanic 
sediments.  Scientists  studying  the  recent 
geological  past  theorize  that  gas-hydrate 
dissociation  during  the  last  glacial  period  (some 
18,000  years  ago)  may  have  been  responsible  for 
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the  rapid  termination  of  the  glacial  episode.  During  the 
glacial  period  the  sea  level  fell  by  more  than  300  feet, 
which  lowered  the  hydrostatic  pressure,  leading  to 
massive  slumping  that  may  have  liberated  large 
volumes  of  methane.  This  may  be  related  to  series  of 
mass  removals  and  glacial  to  non-  glacial  transitions. 
Such  changes  could  have  allowed  methane  to  escape 
into  the  atmosphere  as  a  greenhouse  gas,  possibly 
triggering  climate  change.  As  the  frequency  of 
slumping  and  methane  release  increased,  a  threshold 
was  eventually  reached  where  ice  melting  began, 
leading  to  a  rapid  deglaciation. 

Hydrate  methane  trapped  in  the  sediment  originates 
from  deep  thermogenic  sources  or  is  generated  by 
microbes  in  cycling  of  deposited  organic  matter,  fixed 
through  time  in  deep  sediment.  A  portion  of  the  deep 
sediment  methane  is  transported  by  diffusion  and 
advection  upward  toward  the  overlying  water  column. 
Near  the  sediment-water  column  interface,  a  large 
fraction  of  methane  is  removed  by  microbial  anaerobic 
oxidation  during  which  sulfate  is  reduced.  Since  the 
concentration  of  porewater  sulfate  is  observed  to 
decline  at  or  adjacent  to  the  same  level  in  the  sediment 
where  methane  is  consumed,  it  has  been  proposed  that 
these  two  phenomena  are  linked  through  a  symbiotic 
consortium  of  methanotrophic  archea  and  sulphate- 
reducing  bacteria  (SRB). 

This  shallow  anaerobic  oxidation  of  methane  from 
deep  sediment  provides  an  approach  to  evaluate 
variation  in  the  coastal  methane  vertical  flux  from  deep 
deposits  to  the  shallow  sediment.  Moreover,  from  an 
environmental  and  regulatory  perspective,  it  is 
important  to  understand  the  fundamental  mechanisms 
of  methane  production,  oxidation,  and  transport  in  the 
sediment,  and  the  interplay  between  these  processes 
which  regulates  the  methane  flux  into  the  hydrosphere. 

Currently  seismic  profiling  is  the  primary  method  to 
estimate  the  distribution  of  hydrate  reservoirs  in  the 
seafloor  sediment.  No  other  approach  provides  rapid 
surveys  of  large  regions.  However,  hydrates  are  not 
directly  detected  with  seismic  surveys,  but  the  presence 
is  inferred  by  the  existence  of  a  bottom  simulating 
reflection  (BSR)  with  the  acoustic  transition  through 
the  free  gas  phase.  This  survey  does  not  consisently 
provide  reliable  identification  of  hydrate  beds;  hydrate 
deposits  have  been  incorrectly  identified  or  not 
observed.  To  support  seismic  exploration  of  hydrates, 
geochemical  porewater  profiles  in  shallow  sediments 
are  being  interpreted  to  confirm  deep  methane  deposits. 
Significant  work  is  needed  to  understand  the 


mechanisms  and  interactions  between  transport 
and  methane  sources  and  sinks,  and  to  develop 
models  that  could  be  applied  to  evaluate  the 
methane  distribution  and  assess  environmental 
impacts  of  gas  recovery. 

II.  Objectives 

Research  voyage  Gas  Hydrates  2006  had  9 
primary  objectives  that  are  part  of  the  GNS 
programmes,  National  Petroleum  Resources, 
Offshore  Frontiers,  Gas  Hydrates  (Foundation  of 
Research,  Science,  and  Technology,  FRST, 
Energy  Portfolio),  “How  do  gas  hydrates  weaken 
the  seafloor,  causing  submarine  slides  and 
tsunamis?  ”  (Marsden  Fund,  in  collaboration  with 
NIWA  and  the  University  of  Otago),  and  the 
“DOE  gas  hydrates”  Non  Specific  Output  Fund 
(now  Capability  Fund).  These  objectives  were 
combined  with  NRL  research  focus  on  an 
Advanced  Research  Initiative  (ARI),  Quantify  the 
processes  responsible  for  the  distribution  of  gas 
hydrate  (meters  to  hundreds  of  meters  scale) 
associated  with  representative  seafloor  seeps  and 
an  ONR  sponsored  program  for  the  Marine 
Biogeochemistry  Section,  Naval  Research 
Laboratory  and  Hawaii  Natural  Energy  Institute, 
University  of  Hawaii  for  International 
Collaboration  on  Methane  Hydrate  Research  and 
Development  and  the  Department  of  Energy, 
National  Energy  Technology  Laboratory  National 
Methane  Hydrate  R  &  D.  In  this  broad 
collaboration  the  expedition  also  supported 
research  questions  addressed  by  scientists  from 
the  University  of  Rochester  (New  York,  USA), 
the  GeoForschungsZentrum  Potsdam  (Germany), 
and  the  University  of  Auckland. 

The  voyage  focuses  on  two  study  areas, 
offshore  of  the  Wairarapa  (Study  area  A  or  Ridge 
I)  and  Rock  Garden/Ritchie  Ridge  offshore  of 
Hawke’s  Bay  (Study  area  B).  The  following 
objectives  were  addressed  during  this  research 
expedition. 

Objective  1.  Acquire  high-resolution  seismic 
data  with  coincident  swath  and  3.5  kHz  mapping 
across  a  potential  gas  hydrate  “sweet  spot”  (layer 
with  highly  concentrated  gas  hydrates). 


7 


Objective  2.  Collect  sediment  cores  (piston  or 
gravity)  in  a  transect  across  a  potential  gas  hydrate 
“sweet  spot”  for  pore  water  analysis. 

Objective  3.  Measure  heatflow  in  a  transect  across  a 
potential  gas  hydrate  “sweet  spot”. 

Objective  4.  Collect  microbiological  samples  above 
the  potential  gas  hydrate  “sweet  spot”. 

Objective  5.  Acquire  12  kHz  data  and  measure 
methane  in  the  water  column  using  CTDs  for  detecting 
potential  methane  “plumes”. 

Objective  6.  Retrieve  a  temperature  sensor  deployed 
in  2/2005.  Study  area  B,  Marsden. 

Objective  7.  Collect  sediment  cores  (piston  or 
gravity)  for  reconstruction  of  paleoceanographic 
conditions. 

Objective  8.  Collect  dredge  samples  and/or  sediment 
cores  (piston  or  gravity)  for  investigation  of  rocks  on 
Rock  Garden,  Study  area  B,  Marsden. 

Objective  9.  To  acquire  high-resolution  seismic  data 
with  coincident  swath  and  3.5  kHz  mapping  across 
Rock  Garden  to  study  pockets  of  shallow  gas  beneath 
gas  hydrates. 

This  report  provides  initial  geochemical  and 
heatflow  data.  A  total  project  report  is  available 
through  Dr.  Ingo  Pecher  at  Ingo.Pecher@pet.hw.ac.uk. 

III.  Regional  Description 


freeze-thaw  mechanism  similar  to  frost  heave 
[. Pecher  et  al.,  2005].  This  was  addressed  with 
piston  coring  for  paleoceanographic  temperature 
studies,  additional  seismic  data  acquisition,  and 
coring/dredging  to  survey  sediment  properties. 

Another  part  of  the  expedition  was  to  retrieve 
a  temperature  sensor  on  the  seafloor  that  was 
deployed  in  March  of  2005  in  the  study  area.  The 
combination  of  these  data  with  previously 
collected  data  should  significantly  improve 
scientific  knowledge  of  potential  gas  hydrate 
reservoirs  and  will  allow  a  comprehensive 
analysis  to  test  the  hypothesis  of  “gas  hydrate 
frost  heave”. 


“SS’OO1 
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-39 ’Off 
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Figure  1:  Regional  map  showing  the  two  study  areas. 

Potential  gas  hydrate  desposits,  outlined  by  bottom  Dashed  lines  are  crustal  seismic  data, 

simulating  reflections  in  seismic  data,  appear  to  be 
ubiquitous  on  the  Hikurangi  Margin  [Katz,  1981].  The 
primary  goal  of  this  voyage  was  to  spatially  locate  gas 
hydrate  deposits  as  a  potential  energy  resource.  For  this 
program  available  seismic  data  were  used  to  select 
locations  to  conduct  piston  coring  for  geochemistry  and 
to  measure  heatflow  over  potential  gas  hydrate  “sweet 
spots”,  i.e.,  locations  of  high  gas  hydrate  concentration 
[ Crutchley  et  al ,  2006;  Pecher  and  Henrys ,  2003] 

(Figure  1).  In  addition,  the  voyage  was  intended  to 
contribute  to  data  on  how  gas  hydrates  may  weaken  the 
seafloor  linking  gas  hydrates  with  seafloor  instability, 
submarine  landslides,  and  possibly  Tsunamis  by  a 
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The  broad  plan  of  the  voyage  is  to  undertake 
multi-channel  seismic  surveys,  EM300  multi¬ 
beam,  3.5  kHz,  and  12  kHz  mapping,  sediment 
sampling  by  means  of  drop/piston  coring  and 
rock  dredging,  heatflow  observations  of  the 
seabed,  the  analysis  of  water  masses  and  water 
samples  obtained  by  CTD-drops,  and  a  retrieval 
of  a  temperature  sensor. 

We  are  focusing  on  two  areas.  Ridge  I 
selected  in  the  Study  Area  A  is  located  offshore 
of  the  Wairarapa  (Figure  2),  where  a  suspected 
gas  chimney  has  been  observed  in  seismic  data, 
potentially  indicating  a  deposit  of  highly 
concentrated  gas  hydrates.  Rock  Garden  in 
Study  Area  B  offshore  of  Hawkes  Bay  is  the 
location  at  which  we  suggest  that  gas-hydrate- 
related  seafloor  erosion  occurs  (Figure  3). 

IV.  Cruise  Participants  and  Land-Based 
Collaborators 

Scientists  from  New  Zealand,  United  States, 
Belgium,  Canada  and  Germany,  representing  1 1 
University  and  Government  research  institutions 
contributed  to  this  expedition  (Table  1,  Figure 
4).  Expertise  of  the  science  team  resulted  in  the 
contribution  geochemical,  geophysical, 
geological,  molecular  ecology  and  biological 
data  to  address  ocean  water  column,  sediment, 
and  porewater  research  questions.  Research 
cruise  and  laboratory  participants  are  recognized 
for  the  data  contributions  to  the  research 
questions  outlined  for  this  survey.  This 
diversity  in  expertise,  contributed  through  senior 
scientists,  post  doctoral  researchers,  graduate 
students,  and  technicians  was  necessary  to 
successfully  address  the  research  goals. 


Figure  2:  Study  area  on  Ridge  I  with  piston  coring  and  heatflow 
conducted  on  seismic  Line  1  and  Line  3. 


Figure  3:  Study  Area  B,  Rock  Garden,  dashed  lines  show  most 
relevant  existing  seismic  data,  both  crustal  and  high  resolution. 


Figure  4:  Science  team  participating  in  the  research  cruise  (see 
Table  1). 
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Table  1 :  Cruise  participants,  affiliations  and  project  goals.. 


Pecher 

Dr.  Ingo 

GNS 

Voyage  leader,  seismic, 
multbeam,  3.5  kHz 

Coffin 

Dr.  Richard 

NRL 

Voyage  co-leader,  piston 
core/heat  flow 

Henrys 

Dr.  Stuart 

GNS 

Voyage  co-leader,  seismic, 
multbeam,  3.5  kHz 

Faure 

Dr.  Kevin 

GNS 

Water  column  chemistry 

Greinert 

Dr.  Jens 

GNS/U.  Ghent 
Belgium 

GNS  coordinator  piston 
cores,  water  column 
chemistry,  18  kHz, 
carbonates 

Orpin 

Dr.  Alan 

NIWA 

Paleoceanographic  coring, 
multibeam,  3.5  kHz 

Gerring 

Pete 

NIWA 

Engineer/ technician 

McGill 

Kevin 

NIWA 

Engineer/technician 

Mitchell 

John 

NIWA 

Engineer/ technician 

Wilcox 

Steve 

NIWA 

Engineer/technician 

Gorman 

Dr.  Andrew 

U  of  Otago 

Seismic,  student  advisor 

Crutchley 

Gareth 

U  of  Otago 

Student  watchstander 

Fohrmann 

Miko 

U  of  Otago 

Student  watchstander 

Kilner 

Jeremy 

U  of  Otago 

Student  watchstander 

Senger 

Kim 

U  of  Otago 

Student  watchstander 

Wood 

Dr.  Warren 

NRL 

Heatflow 

Hamdan 

Dr.  Leila 

NRL 

Piston  core/microbiology 

Pohlman 

John 

NRL 

Piston  core/biogeochemistry 

Smith 

Joe 

U  of  Mass  (NRL) 

Piston  core 

Dale 

Celeste 

NRL 

Piston  core 

MacDonald 

Bob 

Biologica  Ltd  (NRL) 

Heatflow/piston  core 

Downer 

Ross 

Milbar  Hydro-Test 
(NRL) 

Piston  core/heatflow, 
engineer 

Bryant 

Layton 

Milbar  Hydro-Test 
(NRL) 

Piston  core/heatflow, 
engineer 

Kurasaki 

Ryan 

U  of  Hawaii 

Microbiology 

Yoza 

Dr.  Brandon 

U  of  Hawaii 

Microbiology 

Tomaru 

Dr.  Hitoshi 

U  of  Rochester 

Iodine  isotopes,  halogens, 
oxygen,  hydrogen 

Kukowski 

Dr.  Nina 

GFZ  Potsdam, 
Germany 

Tectonics,  multibeam 

Schinteie 

Richard 

U  of  Auckland 

Student  watchstander,  pore 
water  analysis,  carbonates 
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V.  Methods 

1 .  Site  Navigation 

The  RV  Tangaroa  (Figure  5)  had  a 
single  propeller,  no  bow  thruster  (during  these 
operations),  and  no  computer  aided  dynamic 
positioning  capability.  The  station-keeping 
accuracy  was  affected  by  the  weather  and 
current  conditions  and  varied  throughout  the 
cruise.  For  accuracy  of  the  sample  station 
selection  it  was  necessary  to  determine  the 
positioning  capability.  This  was  accomplished 
by  monitoring  variation  in  positioning  during 
heatflow  data  acquisition  at  5  stations.  For 
stations  1-3  this  was  approximately  +/-  30m, 
for  station  4  approximately  +/-  10m  and  for 
Station  5  approximately  +/-  50  m.  The  position 
of  the  probe  on  the  seafloor  was  taken  to  be  the 
position  of  the  ship  at  the  time  of  insertion. 
This  variation  in  the  positioning  was 
acceptable  relative  to  the  coring  and  heatflow 
sample  regions  selected  across  the  seismic 
surveys.  Given  the  circumstances,  the  piloting 
of  the  ship  was  extraordinary. 


Figure  5:  RV  Tangaroa  moored  in  Wellington,  New 
Zealand. 

2.  Piston  Coring 

Piston  coring  (PC)  was  conducted  over  the 
mid  starboard  side  of  the  vessel.  PC  was 
conducted  on  a  24  hour  basis  through  a  NIWA 
contract  for  the  ship  costs  (Figure  6).  The  core 
barrel  length  was  6.5  m  and  water  column 
depths  above  the  regions  studied  ranged  from 
670  to  2200  meter.  The  piston  core  trigger  was 


initially  set  at  approximately  3  feet  above 
sediment  penetration.  Cores  retrieved  with  this 
drop  depth  were  too  shallow,  ranging  from  1  - 
3  meters.  An  increase  in  five  feet  for  the 
trigger  drop  depth  increased  the  core  length  to 
3-4  meters  at  most  locations.  When  the 
sediment  structure  did  not  allow  the  core 
penetration  a  2.5  meter  core  barrel  was  used 
and  the  core  weight  was  set  as  a  gravity  core. 


Figure  6:  Piston  coring  system  provided  by  NIWA. 

3.  Pore  Water 

a.  Press  Loading  and  Squeezing  -  Core 
sections  were  retrieved  from  the  piston  core 
barrel  in  two  3  m  sleeves.  10  cm  whole  round 
core  sections,  collected  at  a  vertical  resolution 
of  25  cm,  were  taken  to  the  ship  wet  laboratory 
for  processing  (Figure  7).  Sediment  from  the 
interior  of  the  core  was  placed  inside  the  body 
of  the  press.  The  press  bodies  were  filled  with 
sediment  to  maximum  capacity  with  a 
minimum  air  volume.  A  clean  latex  sheet  was 
placed  on  the  air  inflow  side  of  the  press  body 
to  prevent  free  airflow  through  the  press.  The 
press  bodies  were  capped,  placed  on  the  pore 
water  press  rack  and  pressurized  with  60-80  psi 
air  over  the  latex  sheet.  Pore  water  was 
collected  in  a  60  ml  syringe  on  the  outflow 
line.  Sediment  from  the  press  was  placed  on  a 
sheet  of  ashed  A1  foil  and  frozen. 
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Figure  7:  Porewater  press  lab  setup  with  24  press  ports. 


b.  Pore  Water  Processing  -  The  average 
volume  of  pore  water  collected  for  each  sample 
was  23.4  ml  with  a  range  of  <1  to  38  mis. 

Using  the  side  port  on  a  three-port  stopcock,  1 
ml  was  removed  from  the  60  ml  syringe  with  a 
10  ml  syringe  for  analysis  of  total  sulfides  and 
2  mis  were  removed  for  ion  analysis  (SOT2  and 
Cf)  and  transferred  into  2  ml  screw  top  vials. 

A  volume  of  100  pi  of  0.8  M  CdN03  was 
added  to  the  sample  for  ion  analysis  to 
precipitate  dissolved  sulfides  which  could 
oxidize  to  sulfate.  A  0.2  pm  filter  in  a  holder 
was  placed  on  the  60  ml  syringe  and  the 
remaining  contents  of  the  syringe  were  placed 
into  a  scintillation  vial.  Using  a  1  ml  fixed 
volume  pipet,  2  mis  (exactly)  of  the  sample 
were  transferred  into  a  5  ml  serum  vial  for  DIC 
concentration,  2  mis  were  added  to  a  3  ml 
serum  vial  for  analysis  of  813C-DIC,  813C  - 
DOC  and  DOC  concentration.  The  serum  vials 
were  sealed  with  Teflon  lined  septa.  Samples 
for  sulfide,  ion  and  DIC  concentration  analysis 
were  delivered  to  the  on-board  geochemical 


laboratory  and  processed  immediately.  The 
vials  for  carbon  isotope  and  DOC 
concentration  analysis  were  frozen  at  -20°C. 

Additional  samples  for  analyses  included,  2 
mis  in  5  ml  amber  ampoules  for  DOC 
fluorescence,  2  mis  for  volatile  fatty  acids 
(VFAs)  in  3  ml  glass  screw  top  vials  with 
Teflon  lined  septa,  7  mis  for  nutrients  (NO3, 
NO2,  NH4  and  PO4),  alkalinity,  major  anions 
and  major  cations  (by  ICP)  in  10  ml 
polyethlyene  vials  and  2  mis  for  iodine.  The 
remaining  water  was  archived  in  5  ml  amber 
ampoules. 

4.  Ship  Board  Laboratory  Analysis 

a.  Sulfate  and  Chloride  Concentrations  - 

Sulfate  and  chloride  were  measured  with  a 
Dionex  DX-120  ion  chromatograph  equipped 
with  an  AS-9HC  column  (Figure  8).  Samples 
were  diluted  1:50  (vol/vol)  prior  to  analysis 
and  measured  against  al:50  diluted  IAPSO 
standard  seawater  (28.9  mM  SO4"2,  559  mM 
Cl  ).  Chloride  concentrations  lower  than 
values  measured  in  seawater  is  used  to  interpret 
the  present  of  melted  hydrates  in  the  core.  The 
depth  of  depleted  sulfate  concentrations 
measured  in  porewater  (limits  of  detection  ~0. 1 
mM)  was  used  to  assist  in  the  prediction  of 
hydrate  beds  below  (described  below). 
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Figure  8:  Dionex  IC  used  for  chloride  and  sulfate 
analysis. 

b.  Porewater  Methane  Concentrations  - 
Methane  concentrations  were  determined  from 
3 -ml  sediment  plugs  using  headspace 
techniques  and  were  quantified  against 
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certified  gas  standards  (Scott  Gas, 
Plumbsteadville  PA).  Headspace  analysis  was 
performed  on  board  using  a  GC-FID  Shimadzu 
GC-14A  gas  chromatograph  equipped  with  a 
Hayesep  0.80/100  column  (Figure  9).  Methane 
concentrations  are  presented  in  millimolar  units 
(mM). 


Figure  9:  Shimadzu  GC  14A  used  for  porewater  methane 
analysis. 

c.  Sulfide  concentrations  -  Pore  water  sulfide 
concentrations  were  measured  with 
spectrophotometer  (Figure  10),  using  the  Cline 
method  (Cline,  1969).  Sulfide  concentrations 
are  presented  in  millimolar  units  (mM). 


Figure  10:  ThermoSpectronic  Genesis  20 
spectrophotometer  used  for  porewater  sulfide  analysis. 


d.  Pore  Water  DIC  Concentrations  -  Pore  water 
dissolved  inorganic  carbon  (DIC) 
concentrations  were  measured  using  a  UIC 
coulometer  (Figure  11)  and  standardized 
against  a  certified  reference  material  (CRM, 
Batch  58).  DIC  concentrations  are  presented  in 
millimolar  units  (mM). 


Figure  11:  UIC  Inc.  Coulometric  system  030  used  for 
DIC  analysis. 

5.  Porewater  Data  Interpretation 

Parameters  for  the  porewater  data  include 
methane,  sulfate,  chloride,  sulfide,  and 
dissolved  inorganic  carbon  concentration  and 
stable  carbon  isotope  ratio.  Porewater 
geochemical  data  are  plotted  for  comparison  of 
methane  and  sulfate  concentrations  to  assess 
the  varying  degrees  of  vertical  methane  fluxes 
and  resulting  methane  oxidation  (Borowski  et 
al.,  1999).  Anaerobic  oxidation  of  methane 
occurs  as: 

CH4  +  S042'  -»•  HC03'  +  HS'  +  H20 

in  sediment  depths  where  vertical  flow  of  deep 
sediment  methane  and  sediment  surface  sulfate 
converge.  This  analysis  is  used  to  estimate  the 
presence  of  methane  deep  in  the  sediments  that 
could  be  concentrated  in  hydrate  beds.  DIC  is 
added  in  these  plots  to  support  the 
interpretation  of  methane  oxidation  and  sulfide 
is  incorporated  for  comparison  with  sulfate 
reduction.  Stable  carbon  isotope  analysis  of 
the  DIC  provides  additional  interpretation  of 
methane  oxidation  and  contribution  to  the  DIC. 
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6.  Heat  Flow 

Thermal  data  collected  in  the  upper  few 
meters  of  the  seafloor  using  a  heatflow 
instrument  have  proven  to  be  reliable  for  a 
proxy  for  fluid  flow  and  help  define  the  limits 
of  active  flows  around  methane  seeps  and  mud 
volcanoes  associated  with  methane  seeps  and 
hydrates.  The  heatflow  instrument  used  was  a 
3.5-meter-long  “violin  bow”  or  “Lister-type” 
instrument  (Figure  12).  Eleven  thermistors 
were  arranged  30  centimeters  apart  in  a  1-cm- 
diameter  tube  held  in  tension  parallel  to  a  solid 
steel  strength  member.  There  was  also  a 
temperature  sensor  mounted  on  the  top  of  the 
weight  stand  which  recorded  the  water 
temperature  near  the  sediment-water  interface. 
The  system  measured  both  temperature 
gradient  and  thermal  conductivity  in-situ. 
Sediment  temperatures  were  calculated  from 
the  decay  of  the  frictional  heat  caused  by 
penetration  of  the  instrument  into  the  sediment. 


Figure  12:  Sediment  heatflow  “Lister-type”  system 
deployment. 

Thermal  conductivity  was  determined  from 
the  decay  of  a  calibrated  thermal  pulse  applied 
after  a  preset  period  of  time.  Heat  flow  values 
were  determined  at  each  station  by  computing 
thermal  resistance  values  at  each  thermistor, 

R  =  J  (1/A.)  dz, 


where  A  is  the  thermal  conductivity.  In  a 
situation  of  steady-state  conductivity  the  heat 
flow  is  equal  to  the  slope  of  the  line  on  a 
Bullard  Plot,  a  plot  of  temperature  vs.  thermal 
resistance.  For  each  station,  any  non-linear  data 
that  might  be  attributed  to  bottom  water 
warming  or  cooling,  were  removed  so  as  not  to 
bias  the  statistics.  A  heatflow  value  is 
determined  from  the  slope  of  the  best-fitting 
linear  least-squares  line  through  the  remaining 
data.  All  heatflow  values  were  corrected  for 
instrument  tilt. 

High  resolution  transects  were  done  over 
the  seeps  and  mounds  get  an  accurate  summary 
of  elevated  thermal  signatures.  Stations  were 
typically  set  no  more  than  100  meters  apart  to 
avoid  oversight  on  fluid  flow  associated  with 
distinct  lateral  distribution  of  seafloor  seeps. 
Data  typically  show  clear  anomalies  in 
sediment  temperature  and  heatflow  associated 
with  the  mounds  and  seeps. 

7.  Laboratory  Analysis 

1  'J  1  O 

a.  D1C  8  C  Analysis  -  8  C  is  measured  on 
pore  water  D1C.  D1C  in  porewater  was  released 
as  CO2  into  the  headspace  by  adding  100  pi  of 
10%  HC1  into  sample  serum  vials.  Gases  are 
analyzed  with  a  Finnigan  Delta  S  Isotope  Ratio 
Mass  Spectrometer  with  a  sample  injection 
through  a  Varian  GC.  For  S13C  analysis  the 
standard  is  PeeDee  Belmenite. 

b.  Sediment  porosity  is  calculated  through 
the  core  profiles  with  subcore  sediment 
samples.  Samples  are  placed  in  pre -weighed 
47  mm  Petri  dishes  and  frozen  until  returned  to 
the  laboratory.  Back  in  the  laboratory  samples 
are  thawed,  the  outer  Petri  dish  is  dried  at  room 
temperature  overnight.  Samples  are  then 
placed  in  a  drying  oven  at  50°C  overnight. 
Thoroughly  dried  sediment  is  weighed  again 
and  the  change  is  used  for  the  water  content  in 
the  calculation  of  sediment  porosity  that  is 
incorporated  in  the  calculations  of  vertical 
sulfate  diffusion. 
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8.  Molecular  Biology 

a.  DNA  extraction  and  conditions  for  Length 
Heterogenity-Polymerase  Chain  Reaction  (LH- 
PCR).  Genomic  DNA  was  extracted  from  500 
mg  (wet  weight)  of  sediment  using  Bio  101 
FastDNA  ®  SPIN  kit  for  soil  (Q-BIOgene, 
Inc.),  and  stored  at  -20°C.  Approximately  10 
ng  of  DNA  was  used  as  template  for  LH-PCR 
in  triplicate  reactions.  PCR  amplification  of 
the  first  two  variable  regions  of  the  small 
subunit  (SSU)  rRNA  was  performed  using  a 
fluorescently  labeled  (6-FAM-27F)  forward 
primer  (5’-6-FAM-AGAGTTTGATCM 
TGGCTCAG-3’)  and  a  non-fluorescent  reverse 
primer,  355R  (5’-GCTGCCTCCCGTAGG 
AGT-3’)  (PE  Biosystems).  Positive  (E.  coli 
laboratory  culture)  and  negative  controls 
accompanied  each  reaction.  PCR  mixtures  and 
conditions  are  according  to  Hamdan  et  al.  (in 
prep).  Briefly,  a  final  concentration  of  IX  PCR 
buffer,  2.5  mM  MgS04,  0.20  mM  of  dNTPs, 
0.5U  Taq  Gold  DNA  Polymerase  (Applied 
Biosystems,  Inc.),  0.5  pM  of  primers  and 
DEPC  water  up  to  a  final  volume  of  20  pL  per 
reaction.  A  Gene  Amp  System  9700  (Applied 
Biosystems,  Inc.)  was  programmed  for  an 
initial  denaturation  step  at  95°C  for  15  min. 
followed  by  35  cycles  of  95°C  for  30  sec,  48°C 
for  30  sec.  and  72°C  for  2  min  plus  5sec  per 
cycle  with  a  final  extension  at  72  °C  for  30-45 
min.  LH-PCR  product  was  run  on  a  1% 
agarose  gel  to  visualize  and  quantitate  product. 
The  PCR  Product  was  diluted  based  upon  the 
strength  of  the  product,  mixed  with  internal 
size  standard  (ILS-600,  Promega),  and  HiDi 
formamide  (Applied  Biosystems,  Inc.)  and 
loaded  on  an  SCE  9610  capillary  sequencer 
(Spectrumedix,  LLC).  After  analysis  of  LH- 
PCR  fingerprints  with  GenoSpectrum  software 
(SpectruMedix,  LLC)  a  custom  PERL  script 
was  used  to  extract  the  LH-PCR  data, 
normalize  profiles  with  respect  to  the  total  peak 
area  in  profiles  and  interleave  profiles  being 
analyzed.  Reproducibility  was  tested  by 
comparing  triplicate  PCRs  for  each  sample. 

b.  Cloning  and  sequencing.  New  PCR 
products  with  non-fluorescent  primers  will  be 
cloned  using  a  TOPO  TA  cloning  kit 
(Invitrogen).  Transformed  cells  were  plated  on 


imMedia-AmpBlue  (Invitrogen)  agar  and 
incubated  at  37°C  for  20±4  h.  For  each  clone 
library,  at  least  96  individual  white  colonies 
were  obtained,  and  extracted  DNA  from 
colonies  was  used  as  template  for  PCR  using 
M13  forward  (5'  -  GTAAAACGACGGCCAG 
-  3’)  and  reverse  (5’-  CAGGAAACAGCTAT 
GAC  -  3')  primers  (Invitrogen).  After 
visualizing  on  a  1%  agarose  gel,  PCR  products 
were  purified  with  AMPure  magnetic  beads 
(Agencourt  BioScience  Corp.).  After 
quantitating  purified  products  on  an  agarose 
gel,  sequencing  reactions  were  performed  with 
Big  Dye  Terminator  Cycle  Sequencing  Kit 
v3. 1  (Applied  Biosystemsand  purified  using 
Sephadex  G-50  (Sigma-Aldrich).  Reactions 
were  reconstituted  in  HiDi  formamide  (Applied 
Biosystems)  and  run  on  a  SEC9610  capillary 
sequencer  (SpectruMedix  LLC).  After 
analyzing  raw  data  with  BaseSpectrum 
software  (SpectruMedix  LLC),  sequences  were 
imported  in  Sequencher  (Gene  Code  Corp.)  for 
manual  base  calling. 

Clone  sequence  data  were  used  in 
MEGABLAST  searches  against  GenBank,  and 
the  Ribosomal  Database  Project  (RDP)  v8.1. 
Clone  libraries  were  analyzed  using  classic 
rarefaction  analysis  to  determine  the  estimated 
diversity  of  the  phylotypes  in  the  community. 
Clone  sequences  were  aligned  with  matching 
sequences  from  GenBank  or  RDP  v8.1  using 
Clustal  X  and  phylogenetic  trees  were 
constructed  using  PAUP. 

c.  Nucleotide  sequence  accession  numbers. 
Partial  environmental  16S  rDNA  sequences 
were  deposited  in  the  GenBank  database  and 
be  made  available  to  the  public. 

VI.  Preliminary  Geochemical  Results 

The  following  presentation  of  geochemical 
and  heat  flow  data  will  be  finished  with 
laboratory  sample  and  data  analyses.  These 
data  are  presented  for  overview  of  the 
preliminary  results  and  to  provide  focus  for  the 
key  issues  in  data  interpretation.  These  data 
will  be  compared  to  recent  and  current  field 
activity  that  this  research  team  is  involved  with 
on  the  mid  Chilean  coast,  Cascadia  Margin  off 
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the  western  coast  of  Canada,  and  Texas- 
Louisiana  Shelf  in  the  Gulf  of  Mexico. 

A  core  log  summary  is  presented  in  Table 
2.  These  cores  were  obtained  at  a  nearshore 
region  to  the  northeast  of  Wellington,  NZ, 
across  Ridge  I  and  over  the  Rock  Garden,  on 
the  Hikurangi  Margin.  Data  presented  in  this 


report  are  the  first  geochemical  field  analysis  to 
survey  methane  hydrate  deposits  along  the 
coast  of  New  Zealand.  These  data  will  be 
available  for  future  methane  hydrate 
exploration  in  this  region. 


Table  2:  Piston  coring  log  for  samples  taken  at  the  neashore  control  site,  Ridge  I  and  Rock  Garden. 


Cruise  log  ID 

NRL  log  ID 

Seismic 

line 

wc 

depth 

(m) 

Core 

depth 

(cm) 

Lat-long 
(40  meters  up) 
(Core  in) 

Pull 

out 

(lbs) 

Comment 

Lab 

Nearshore 

TAN0607-2 

NZ  Core  1 

S40°93.687 

E176°72.492 

NIWA 

TAN0607-3 

NZ  Core  2 

687 

350 

S40°93.687 

E176°72.492 

shallow  test 
and  control 

core 

NRL 

Ridge  I 

TAN0607-17 

NZ  Core  3 

4610 

1998 

332 

S40°50.143 

E177°22.242 

9.4 

sleeve  cracked 
the  full  length 

NRL 

S40°50.143 

E177°22.242 

TAN0607-18 

NZ  Core  4 

4510 

480 

S40°50.498 

E177°22.518 

9.2 

NRL 

S40°50.501 

E177°22.509 

TAN0607-19 

NZ  Core  5 

4098 

2128 

— 

S40°5 1.408 
E177°25.993 

8.0 

core  did  not 
trigger 

NRL 

S40°51.414 

E177°25.994 

TAN0607-20 

NZ  Core  5 

4098 

2124 

285 

S40°5 1.406 
E177°24.990 

10.3 

NRL 

S40°5 1.404 
E177°25.002 

TAN0607-22 

NZ  Core  6 

4279 

1912 

— 

S40°5 1.009 
E177°24.471 

10 

bent  pipe  no 
core  retrieved 

NRL 

S40°51.017 

E177°24.468 

TAN0607-23 

NZ  Core  7 

4230 

2020 

275 

S40°5 1.053 
E177°24.669 

10.2 

sulfide  odor 

NRL 

S40°5 1.030 
E177°24.691 

TAN0607-24 

NZ  Core  8 

4255 

1993 

250 

S40°51.116 

E177°24.872 

9 

NRL 

S40°51.118 

E177°24.884 

TAN0607-25 

NZ  Core  9 

4400 

1993 

S40°50.751 

E177°23.449 

9.2 

Light  sulfide 

NIWA 

TAN0607-26 

NZ  Core  10 

4175 

2068 

98 

S40°51.135 

E177°25.336 

10.2 

no  pipe  bend, 
stiff  sediment 

NRL 

S40°5 1.237 
E177°25.340 

TAN0607-27 

NZ  Core  11 

4400 

1981 

250 

S40°50.735 

9.2 

lost  trigger 

NRL 

16 


E177°23.445 

core 

S40°50.737 

E177°23.466 

TAN0607-28 

NZ  Core  12 

4371 

1945 

128 

S40°50.804 

E177°22.688 

9.6 

trigger  cable 
extended 

NRL 

S40°50.799 

E177°23.705 

Bent  core 

TAN0607-45 

NZ  Core  13 
(20) 

4175 

2089 

240 

S40°51.2352 

E177°25.3319 

10.0 

core  logging 
changed  by  I. 
Pecher  in  the 

NZ  notation. 
Core  #’s  in 
their  log  are 
NRL  numbers 
+  7. 

NRL 

S40°5 1.2346 
E177°25.3321 

TAN0607-46 

NZ  Core  14 
(21) 

4200 

2031 

450 

S40°51.1700 

E177°25.1399 

9.8 

statement  from 
first  mate 
navigation 
resolution  is 

50  m 

NRL 

S40°51.1718 

E177°25.1459 

TAN0607-47 

NZ  Core  15 
(22) 

4135 

2134 

424 

S40°51.3217 

E177°25.6815 

9.6 

NRL 

S40°5 1.3225 
E177°25.6837 

TAN0607-48 

NZ  Core  16 
(23) 

4400 

1979 

— 

S40°50.7377 

E177°23.4594 

8.9 

core  cap  was 
lost  -  no  core 

NIWA 

S40°50.7396 

E177°23.4554 

TAN0607-49 

NZ  Core  17 
(24) 

4385 

1968 

220 

S40°50.7605 

E177°23.5724 

10.3 

Bent  core  pipe 

NRL 

S40°50.7710 

E177°23.5772 

TAN0607-50 

NZ  Core  18 
(25) 

4330 

1912 

200 

S40°50.8910 

E177°24.0255 

8.3 

gravity  core 
used  on  top  of 
the  mound 

NRL 

S40°50.8970 

E177°24.0267 

TAN0607-51 

NZ  Core  19 
(26) 

4330 

1911 

NIWA 

S40°50.8777 

E177°24.0300 

12.0 

Short  ~5M 
piston  core 
used 

NIWA 

S40°50.8889 

E177°24.0218 

Bent  core  pipe 

TAN0607-73 

NZ  Core  28N 

4155 

2119 

310 

S40°51.2817 

E177°25.5062 

9.8 

Last  core  on 
this  line 

NRL 

S40°5 1.2799 
E177°25.51 19 

Trigger  was 
set  back  to 
previous 
length 

TAN0607-74 

NZ  Core  29N 

1319 

1962 

350 

S40°48.847 

E177°24.351 

9.2 

First  core  2  km 
north 

NRL 

S40°48.8461 

E177°24.3447 

TAN0607-75 

NZ  Core  30N 

1405 

1952 

381 

S40°48.930 

E177°24.7130 

9.5 

NRL 

S40°48.9400 

E177°24.7161 

17 


TAN0607-76 

NZ  Core  3  IN 

1536 

1987 

S40°49.0742 

E177°25.2670 

9.6 

NRL 

S40°49.0704 

E177°25.2765 

TAN0607-77 

NZ  Core  32N 

1777 

2086 

400 

S40°49.3512 

E177°26.2712 

8.8 

NRL 

S40°49.3561 

E177°26.2804 

Rock 

Garden 

TAN0607-62 

NZ  Core  22 

891 

350 

S40°03.66 

E178°08.47 

NRL 

TAN0607-63 

NZ  Core  23 

892 

60 

S40°043.51 

E178°87.03 

Very  thick 
clay 

NRL 

TAN0607-64 

NZ  Core  24 

2523 

350 

S40°06.6466 

E178°21.2087 

10.7 

NRL 

S40°06.6457 

E178°21.2079 

TAN0607-65 

NZ  Core  25 

2693 

NIWA 

S40°05.5965 

E178°18.8682 

9.3 

Short  core 
redo 

NIWA 

S40°05.6966 

E178°18.8679 

The  trigger 
chain  was 
shortened 

TAN0607-66 

NZ  Core  26 

2693 

NIWA 

S40°05.5736 

E178°18.8488 

No  core  bent 
barrel 

NIWA 

S40°05.5736 

E178°18.8436 

A.  Core  Testing  Region 

Preliminary  instrument  and  equipment 
testing  was  conducted  at  a  shallow  region  (670 
m)  outside  Wellington  S40°48.02,  W176°54.5 
(Figure  13).  Porewater  sulfate,  methane  sulfide 
and  DIC  concentrations  from  a  core  in  this 
region  provided  a  control  data  base  away  from 
methane  hydrate  deposits  (Figure  14).  Seismic 
surveys  indicated  no  bottom  simulating 
reflector  and  hydrates  were  not  present  in  the 
sediment.  The  water  column  depth  in  this 
region  was  687  m.  The  piston  core  retrieved 
was  350  cm  long.  Sulfate  profiles  showed  a 
shallow  slope  through  the  core  that 
corresponded  with  near  zero  methane 
concentrations.  The  sulfide  and  DIC 
concentrations  show  conservative,  slight 
increases  through  the  core  concentrations. 
These  control  site  data  will  be  compared  with 
SMIs  and  related  geochemical  data  over  Ridge 
I  and  Rock  Garden. 


Line  05CM-3G,  warm-up  core 


Figure  13:  Seismic  line  for  the  control  piston  coring 
(indicated  by  the  arrow)  and  heatflow  analysis.  CDP 
spacing  was  set  at  12.5  m. 
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Figure  14:  Porewater  sulfate,  DIC,  sulfide  and  methane 
concentrations  from  a  piston  core  taken  at  a  non  active 
control  region. 

B.  Ridge  I 

1.  Geochemistry/Piston  Coring 

Piston  coring  was  conducted  along  the 
seismic  profile  2  on  Ridge  I  (Figure  15). 
Seismic,  geochemical  analysis  of  porewaters 
and  heatflow  are  combined  for  a  survey  of 
potential  deep  hydrate  deposits  in  this  study. 
Coring  sites  selected  were  planned  to  assess 
interpretation  of  the  seismic  profiles  for  strong 
BSR  and  vertical  fluid  migration  to  the 
sediment  surface  and  high  gas  wipe  out 
regions.  Thorough  seismic  interpretation  is 
being  prepared  by  I.  Pecher  (GNS). 


Una  0SCM-3a,  coring  and  haatttow  Iraraed 
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Figure  15:  Ridge  I  piston  core  sites  presented  along  the 
seismic  CDP  number.  CDP  spacing  was  set  at  12.5  m. 

Figures  presented  here  are  intended  for  a 
preliminary  overview  of  the  data  acquired 
during  this  expedition.  The  total  data  set  is 


presented  in  Appendix  1 .  Plots  are  included  in 
Figure  16  for  an  outline  the  vertical  methane 
diffusion  in  the  sediment,  in  terms  of  methane 
and  sulfate  profiles,  on  the  landward  side, 
seaward  side  and  top  of  the  mound.  Overall  the 
sulfate  profiles  from  these  piston  cores  were 
non-conservative.  Non-conservative  sulfate 
profiles  in  porewaters  can  result  from  several 
different  factors;  including  variation  in  vertical 
diffusion  caused  by  recent  deposits  of  modem 
sediments,  variation  in  vertical  sediment 
properties,  horizontal  porewater  advection,  and 
bioturbation.  Observation  of  these  cores  found 
worm  channels  in  the  surface  sediment  that 
suggests  bioturbation  causes  vertical  methane 
profiles  in  the  shallow  core  regions.  Vertical 
sulfate  diffusion  was  calculated  on  the  linear 
points  of  the  sulfate  profile  below  segments 
where  bioturbation  appeared  to  mix  seawater 
sulfate  through  the  core. 

A  comparison  of  the  vertical  variation  in 
methane  and  sulfate  (Figure  16)  results  in 
inverse  concentration  gradients,  with  declines 
in  sulfate  from  reduction  to  sulfide  coinciding 
with  increases  in  the  DIC  upon  anaerobic 
oxidation  of  methane.  Core  3  was  taken  on  the 
landward  side,  well  off  the  mound,  in  a  region 
with  a  level  BSR  observed  in  the  seismic 
profile  (Figure  16).  The  porewater  sulfate 
profile  in  this  core  was  similar  to  the  data  from 
the  control  site  with  a  shallow  slope.  Methane 
concentrations  in  the  porewater  from  this 
location  were  at  the  limits  of  detection.  Similar 
sulfate  and  methane  profiles  were  measured  in 
Core  1 8  on  top  of  the  mound  located  in  the 
Ridge  I  survey  area.  These  data  indicate  low 
vertical  porewater  and  methane  diffusion, 
perhaps  impeded  by  carbonate  formation  on 
the  mound.  Rapid  vertical  methane  diffusion 
was  observed  in  the  sulfate  and  methane 
profiles  at  the  coring  locations  on  the  landward 
and  seaward  edges  of  the  mound.  Cores 
obtained  near  the  mound  were  consistently 
shallower  and  methane  concentrations  through 
the  SMI  were  depleted  to  the  limits  of 
detection.  The  sulfate  and  DIC  concentrations 
indicate  that  this  was  a  region  in  which 
anaerobic  oxidation  of  methane  was  active. 
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Sulfate  and  Methane  Profiles 
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Figure  16:  Ridge  I  piston  core  sulfate  and  methane  porewater  data  presented  along  the  seismic  CDP  number.  CDP  spacing 
was  set  at  12.5  m. 


A  summary  of  DIC  profiles  in  sediment 
porewaters  are  presented  in  Figure  17.  The 
concentration  of  DIC  follows  the  sulfate 
profiles  with  the  variation  in  the  vertical 
methane  diffusion  across  the  seismic  line. 
Shallow  DIC  concentration  gradients  were 
measure  on  top  of  the  mound  and  at  the 
landward  site  off  the  mound.  Gradients  on  the 
landward  and  seaside  slopes  of  the  mound  were 
significantly  steeper. 

Sulfide  concentrations  in  porewaters  follow 
the  variations  in  DIC  resulting  from  methane 
oxidation  during  sulfate  reduction  (Figure  18). 
Sulfide  concentrations  were  below  the  limits  of 


detection  on  the  mound  and  at  the  landward 
site.  An  interesting  comparison  of  the  sulfide 
profiles  is  the  large  difference  in  the 
concentration  gradients  between  the  landward 
and  seaward  slopes  (Figure  18).  While  there 
was  not  a  large  difference  in  sulfate  profiles 
between  these  core  sites,  the  higher  sulfide  end 
product  in  the  sulfate  reduction  suggests  a 
greater  anaerobic  oxidation  of  methane  at  the 
landward  coring  locations.  This  observation 
supports  the  heatflow  data  presented  below 
with  active  fluid  advection  on  the  landward 
side  of  the  mound  and  diffusion  on  the  seaward 
point  (presented  below). 
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DIC  and  Sulfate  Profiles 
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Figure  17:  Porewater  DIC  and  sulfate  concentrations  for  samples  representative  of  variations  during  coring  on  Ridge  I. 
CDP  spacing  12.5  m. 
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Figure  18:  Porewater  sulfide  concentrations  at  piston  core  sites  through  Ridge  I.  Total  data  are  presented  and  a  comparison 
of  cores  1 3  and  1 7  is  included  for  additional  data  interpretation. 
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Figure  19:  SMI  calculated  for  porewaters  from  piston  cores  at 
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Figure  20:  Porewater  geochemical  profiles  selected  from  locations  on  the  seaward  (Core  32)  and  landward  (Core  29)  side 
of  the  Ridge  I  mound. 


SMI  measured  across  the  Ridge  I  seismic 
profile  is  presented  in  Figure  19.  The  deepest 
measured  SMI  was  at  the  nearshore  control 
region  at  39.5  m  (Figure  13).  On  the  landward 
side  of  the  Ridge  I  transect  the  deepest  SMI  was 
12.9  m.  Moving  toward  the  slope  of  the  mound 
the  SMI,  resulting  from  a  steep  sulfate  gradient, 


was  1.84  m  in  Core  17.  A  similar  SMI,  2.11  m, 
was  measured  on  the  seaward  side  of  the  mound 
at  Core  7.  Geochemical  porewater  profiles  did 
not  provide  a  sulfate  gradient  on  top  of  the 
mound  (C12). 

Another  transect  was  taken  on  a  parallel 
seismic  survey  along  Ridge  I  (Figure  2).  Time 
allowed  for  4  cores  to  be  taken  in  this  region. 
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Rock  Garden  Coring  Locations 


The  slopes  of  the  sulfate,  H2S,  CH4  and  DIC 
profiles  were  considerably  steeper  on  the 
landward  side  of  this  coring  transect  (Figure  20). 
These  data  were  similar  in  the  spatial  variation  in 
comparison  to  the  primary  Ridge  I  data  base. 

The  SMI  measured  through  this  location  ranged 
from  3.34  to  9.45  m  below  the  seafloor  (Figure 
21).  Similar  to  the  previous  coring  line  the  SMI 
was  shallow  on  the  landward  side  and  deeper 
SMIs  were  measured  on  the  seaward  side  of  the 
mound.  The  pore  water  H2S  concentrations  were 
also  greater  on  the  landward  side,  as  observed 
along  the  primary  geochemistry  transect. 

Pore  water  Sulfate  Methane  Interface 


Figure  21 :  Overview  of  coring  locations  and  SMI 
measured  close  and  parallel  to  the  primary  coring  on  the 
Ridge  I  seismic  profile. 

C.  Rock  Garden 

Preliminary  coring  was  conducted  on  the 
landward  and  seaward  edges  of  a  mound  located 
in  Rock  Garden  (Figure  22).  Porewater  sulfate, 
methane,  DIC  and  sulfide  concentration  profiles 
were  similar  to  values  measured  along  the  Ridge 
(Appendix  I).  The  SMI  estimated  for  Cores  22 
and  24  was  341  cm  and  398  cm,  respectively 
(Figure  23).  Future  fieldwork  will  provide  time 
for  further  investigation  of  this  area. 


Figure  22:  Seismic  profiles  and  piston  core  locations  on 
Rock  Garden 
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Figure  23:  Porewater  geochemical  profiles  for  cores  taken 
on  the  seimic  line  across  Rock  Garden. 
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VII.  Preliminary  Heatflow  Results 

A.  Background 

The  thermometry  program  on  NZ06  was 
designed  to  transect  a  fold  of  accreted  sediments 
on  the  Hikurangi  Margin  (Figures  1  &  2). 
Reflection  seismic  profiles  of  this  particular  fold 
exhibit  a  bottom  simulating  reflector  (BSR) 
indicative  of  methane  hydrate,  and  show 
relatively  undeformed  slope  fill  sediments  on 
both  the  landward  and  seaward  side  of  the  NE- 
SW  trending  fold,  which  is  largely  opaque 
seismically.  Numerous  faults  and  gas  wipeouts 
exhibited  in  the  seismic  profiles,  together  with 
the  BSR,  suggest  significant  fluid  flux  at  this 
location  on  the  Margin,  and  it  is  this  presumed 
flux  (and  the  resulting  thermal  anomaly)  that  was 
the  target  of  investigation  along  the  two,  parallel 
thermometry  dip  transects  across  the  fold. 

B.  Operations 

The  thermometry  data  were  acquired  with  a 
violin-bow  type  instrument  rented  from  the 
Canadian  Geological  Survey  (PGC)  and  operated 
under  contract  by  Robert  MacDonald.  The  probe 
contains  eleven  independent  seafloor  thermistors, 
and  one  water  thermistor  at  the  top  of  the  case 
yielding  twelve  independent  temperature 
measurements  from  0  to  3.3  meters  below  the 
seafloor. 

The  data  were  acquired  in  five  deployments 
(stations)  varying  from  two  to  eleven 
penetrations  per  station.  A  total  of  forty 
penetrations  were  acquired,  of  which  two  (HF38, 
and  HF39)  were  discarded  due  to  excessive  tilt. 
For  each  penetration  the  thermal  probe  was 
inserted  into  the  sediments  for  7  minutes  after 
which  a  heat  pulse  was  triggered.  Extraction 
occurred  five  to  seven  minutes  after  the  heat 
pulse  to  allow  sufficient  decay  to  calculate 
thermal  conductivity  and  subsequently  the  total 
heat  flux.  The  final  calculations  for  thermal 
conductivity  and  heat  flow  will  be  done  on  shore, 
but  preliminary  calculations  of  the  steady-state 
temperature  were  performed  via  a  method 
described  below.  It  is  these  preliminary 
temperatures  that  are  displayed  in  the  figures  in 
this  report  (Figures  24  &  25).  Each  penetration  is 
named  uniquely,  with  the  designation  HF 


followed  by  a  sequential  number,  such  that  the 
penetrations  are  labeled  HF01-HF40. 


NZ06  Stations  1-5  (preliminary) 


Figure  24:  Heatflow  profiles  measured  at  the  control  site, 
seaward  and  landward  of  the  mound. 
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black  -  penetration  1 
red  -  penetration  2 
green  -  penetration  3 
blue  -  penetration  4 
purple  -  penetration  5 
yellow  -  penetration  6 
orange  -  penetration  7 
light  gray  -  penetration  8 
medium  gray  -  penetration  9 
dark  gray  -  penetration  10 


Figure  25:  Heatflow  profiles  measured  across  the  seismic 
profile  located  2  km  north  of  the  primary  seismic  line. 

C.  Acquisition 

Station  1  was  used  to  test  ship  board 
deployment,  and  instrument  operation,  as  well  as 
to  serve  as  a  background  site  to  determine  the  no- 
flux  thermal  gradient  in  this  region. 

Temperatures  are  rather  high  in  this  region 
(about  7  Deg.  C)  due  to  the  shallow  location  (721 
meters). 

Station  2  (penetrations  1-5,  HF11-HF15), 
Station  3  (penetrations  1-10,  HF16-HF25),  and 
Station  5  (penetrations  1-4,  HF37-HF40)  were  all 
acquired  along  TAN0607  Seismic  Line  2  (See 
Figure  26). 

Station  4  (penetrations  1-11,  HF26-HF36) 
was  acquired  along  TAN0607  Seismic  Line  3 
(See  Figure  27).  Aside  from  the  two  discarded 
profiles  the  thermal  data  are  remarkably 
consistent  between  neighboring  sites,  suggesting 
that;  1)  the  profiles  are  adequately  spatially 
sampled,  i.e.  the  major  character  of  the  thermal 
regime  is  likely  to  have  been  captured  in  our 


measurements,  and  2)  the  data  were  properly 
acquired,  with  small  random  errors. 


TAN0607  Thermometry  Line  2  (US.  NRL) 


2.39-2.42 

2.42-2,45 
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TAN 0607  Sesimic  Line  2  (New  Zealand  GNS) 


Figure  26:  Heatflow  data  interpretation  across  TAN0607 
seismic  line  2. 


D.  Preliminary  Processing 

The  initial  processing  was  performed  using 
the  software  supplied  with  the  thermal  probe, 
such  that  data  were  reduced  into  *.pen  files.  The 
software  used  to  analyze  the  *.pen  files  for 
steady  state  temperature  and  thermal 
conductivity  was  not  provided  so  preliminary 
estimates  were  based  on  fitting  the  decay  curves 
in  each  of  the  *.pen  files  with  curves  of  the  form; 

(T-T0)/(  Th  -  T0)  =  >/2  (erf[(r+w)/ (2sqrt(kt))]  - 
erf[(r- w)/ ( 2  sqrt(  kt) )] }  { erf[(  w)/(  2sqrt(kt ))] } 

where  T  is  the  temperature  at  any  given  time  (t) 
and  radius  (r)  from  a  cylindrical  rod  of  diameter 
(w)  in  a  homogeneous  medium  with  a  thermal 
conductivity  (k).  Th  and  T0  are  the  high  and  final, 
steady  state  temperature  respectively.  A  fit  of 
this  equation  to  the  observed  curve  was  made  by 
minimizing  the  difference  between  the  observed 
thermal  decay  and  the  curve  above  computed  for 
various  selections  of  k,  Th  and  T0  (assuming 
w=1.0  cm,  and  r=0).  The  selections  were  guided 
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by  a  simulated  annealing  process  in  which  many 
tens  of  thousands  of  parameters  could  be  tried 
and  only  the  best  selected  as  the  best  fit,  (Figure 
28).  The  resulting  steady  state  temperature  T0 
was  used  in  the  temperature  plots  in  this  report. 
All  the  data  files,  as  well  as  the  programs  used  to 
process  them  are  included  on  the  disk  with  this 
report.  Shore  based  processing  with  standard 
thermal  processing  software  is  anticipated. 


TAN0607  Thermometry  Line  3  (U  S.  NRL) 


TAN0607  Sesimic  Line  3  (New  Zealand  GNS) 


Figure  27:  Heatflow  data  interpretation  across  TAN0607 
seismic  line  3 
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Figure  28:  Modeling  and  data  collection  comparison  of 
heatflow  temperature  pulses 


E.  Preliminary  Interpretation 

Figures  26  and  27  clearly  show  a  major 
thermal  perturbation  associated  with  the 
accretionary  margin  fold.  However,  unlike  a 
single  seafloor  seep,  there  does  not  appear  to  be 
significant  flux  associated  with  the  fold  itself. 
Most  of  the  thermal  profiles  exhibit  linear 
character  (Figures  24  &  25)  suggesting  minimal 
if  any  fluid  advection.  The  exception  is  the 
profiles  on  the  landward  side  of  the  ridge  on  Line 
2,  (Figure  24)  whose  concave  downward 
character  is  indicative  of  fluid  advection. 

The  second  obvious  feature  of  the  thermal 
profiles  is  the  presence  of  a  laterally  extensive 
cool  zone  in  the  regions  dominated  by 
conduction  (all  of  Line  3  and  the  seaward  portion 
of  Line  2).  The  cause  for  this  is  unknown,  but  it 
is  unlikely  an  artifact  of  processing,  as  it  is  not 
present  in  all  the  data.  Quantitative  estimates  of 
fluid  flux  await  final  processing  of  the  data  to 
obtain  values  of  thermal  conductivity  and  heat 
flux. 

VIII.  Preliminary  Molecular  Biological 
Results 

NRL  cores  4,  7,  28  and  32  from  Ridge  I  and 
core  22  from  Rock  Garden  were  selected  for 
molecular  biological  analysis  based  upon  the 
sediment  and  porewater  geochemical 
characteristics.  The  SMI  depth  was  specifically 
targeted  for  molecular  biological  analysis.  This 
limited  core  selection  to  locations  where  the  SMI 
was  shallow  enough  such  that  the  core  barrel  was 
able  to  penetrate  the  SMI.  SMI  depths  for  Ridge 
I  cores  are  provided  in  Figure  19.  The  SMI 
depth  for  core  22  was  3.41m  which  was  nearly  at 
the  bottom  of  the  core.  With  respect  to  methane 
concentrations,  cores  4,  7  and  28  were  active  flux 
cores  whereas  cores  22  and  32  represented 
control  cores  where  methane  flux  was 
significantly  reduced. 

DNA  was  purified  and  extracted  from  68 
samples  upon  return  to  the  NRL  laboratory. 
Extracted  DNA  was  sufficiently  low  in 
concentration  such  that  it  could  not  be  visualized 
on  an  agarose  gel  prior  to  amplificantion.  Initial 
attempts  to  amplify  genomic  DNA  failed  to  yield 
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PCR  product.  This  is  due  to  the  presence  of 
unidenfied  inhibitory  substances  which  could  not 
be  separated  from  purified  DNA.  The  problem 
was  overcome  by  diluting  purified  DNA  1:100 
with  DEPC  water  prior  to  introducing  it  into 
PCR  reaction  mixtures.  Three  replicate  PCR 
reactions  were  performed  on  each  sample 
yielding  variable  results.  The  amplified  products 
are  depicted  in  Figure  29. 


Figure  29:  PCR  product  of  68  samples  of  purified  genomic 
DNA  visualized  on  1%  agarose  gel  by  staining  with 
ethidium  bromide  and  viewed  under  UV.  The  intensity  of 
bands  represents  the  strength  of  the  produce  in  each 
sample.  Positive  and  negative  controls  (far  right)  as  well 
as  a  DNA  sizing  ladders  (far  left)  accompany  each  set  of 
reactions. 


Samples  which  had  sufficient  PCR  product  to 
visualize  in  three  replicate  samples  were 
analyzed  by  LH-PCR.  The  result  of  this  analysis 
was  peaks  of  variable  amplicon  sizes  in  an 
electropherogram,  each  representing  an 
individual  phylotype.  In  addition,  the  intensity 
of  the  peak  (peak  area)  approximates  the 
abundance  of  specific  amplicons. 

Comparison  of  previous  data  obtained  with 
this  method  from  other  locations  (mid-Chilean 
margin),  bacterial  diversity  in  Hikurangi  Margin 
sediments  was  high  with  on  average  1 8 
amplicons  per  sample  compared  to  an  average  14 
amplicons  per  sample  from  the  mid-Chilean 
margin.  In  previous  works,  amplicon  373.6  was 
observed  in  nearly  all  samples.  This  observation 
was  repeated  in  the  study  (Figure  30).  Although 
this  amplicon  corresponded  to  several  different 
taxa  (Hamdan  et  al.  in  prep)  it  is  likely  that  some 
will  match  taxa  observed  during  this  study  when 
sequencing  analysis  is  performed  on  selected 
samples. 


The  intensity  of  bands  which  represents  the 
PCR  reaction  efficiency  was  strongest  in  samples 
from  the  SMI.  In  particular  sample  30  depicted 
in  Figure  29  represents  the  sample  collected  at 
the  SMI  in  NRL  core  4.  Likewise,  sample  59 
which  had  the  strongest  overall  band  intensity 
was  collected  at  the  SMI  from  NRL  core  7. 
Finally,  sample  243  from  one  of  the  lower  flux 
cores  had  a  strong  band  signature  for  the  sample 
collected  at  3.10m,  near  the  SMI. 


Figure  30:  LH-PCR  electropherogram  of  PCR  product 
from  NRL  core  7  SMI  sample  (Sample  59).  The  intensity 
of  peaks  represents  the  abundance  of  each  amplicon  and  its 
corresponding  DNA  in  each  sample.  Amplicon  357.  6 
(358.0)  was  the  largest  peak  in  three  replicate  PCR 
reactions  indicating  that  it  dominates  the  bacterial 
community  in  this  sample. 
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During  a  previous  study,  we  identified  a  single 
amplicon  length,  358.0  which  dominated  SMI 
community  profiles  in  several  cores  over  a 
methane  seep  and  one  over  a  hydrate  mound. 

This  amplicon  was  seldom  detected  in  the  data 
set  for  this  study.  However  a  large  peak 
corresponding  to  this  amplicon  was  detected  in 
sample  30  from  the  SMI  in  NRL  core  4. 

Likewise,  this  amplicon  was  the  largest  peak 
detected  in  sample  59  which  corresponds  to  the 
SMI  sample  from  NRL  core  7  (Figure  30).  This 
peak  accounts  for  approximately  47%  of  peak 
area  for  this  core  indicating  that  the  phylotype 
corresponding  to  this  single  amplicon  nearly 
accounts  for  the  majority  of  bacterial  diversity. 
Based  upon  our  previous  investigation  it  is  likely 
that  the  phylotype  associated  with  this  amplicon 
corresponds  to  a  relative  of  the  Desulfosacccinia 
variabilis  subgroup  of  sulfate  reducing  bacteria 
which  have  been  documented  to  be  key  players 
in  the  anaerobic  oxidation  of  methane  at  the  SMI 
in  numerous  locations  around  the  world 
(Hamdan  et  al.  in  prep  and  references  therein). 
These  data  indicate  that  the  LH-PCR  analysis  is  a 
highly  reliable  tool  for  pinpointing  the  exact 
location  of  the  SMI  and  potentially  the  location 
where  anaerobic  oxidation  of  methane  is 
primarily  occurring. 

The  phylogenetic  analysis  will  proceed  with 
the  cloning  and  sequencing  of  SMI  samples  as 
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APPENDIX  I:  Porewater  data  log. 


Hikurangi  Margin,  TAN0706/NZ0606,  June- July  2006 
_ Pore  Water  and  Gas  Data _ 

Pore  Gas 

Water  Sample 

Sed  Sediment 


Core 

Depth 

(cm) 

Chloride 

(mM) 

Sulfate 

(mM) 

Sulfide 

(mM) 

DIC 

(mM) 

Depth 

(cm) 

ch4 

(mM) 

c2h6 

(nM) 

c,/c2 

(vol) 

SMI 

(cm) 

2 

5 

10 

0.0003 

bd 

4987 

2 

45 

543.2 

27.59 

0.00 

2.20 

50 

0.0006 

bd 

2 

70 

541.0 

27.29 

0.00 

3.41 

75 

0.0006 

bd 

2 

95 

540.8 

27.45 

0.00 

3.52 

100 

0.0007 

bd 

2 

120 

541.2 

27.19 

0.00 

3.70 

125 

0.0008 

bd 

2 

145 

538.4 

27.13 

0.00 

3.80 

150 

0.0008 

bd 

2 

170 

543.9 

27.12 

0.00 

3.97 

175 

0.0009 

bd 

2 

195 

545.2 

27.08 

0.00 

3.99 

200 

0.0009 

bd 

2 

220 

542.7 

27.04 

0.01 

4.10 

225 

0.0011 

bd 

2 

245 

560.8 

27.67 

0.05 

4.14 

250 

0.0012 

bd 

2 

270 

540.2 

26.58 

0.06 

4.18 

275 

0.0010 

bd 

2 

295 

541.4 

26.11 

0.00 

4.28 

300 

0.0012 

bd 

2 

320 

541.2 

26.15 

0.30 

4.40 

325 

0.0014 

bd 

2 

345 

548.3 

25.98 

0.06 

4.56 

350 

0.0014 

bd 

3 

28 

544.3 

27.62 

0.01 

2.95 

33 

0.0001 

bd 

1287 

3 

57 

544.7 

27.41 

0.01 

3.34 

63 

0.0001 

bd 

3 

87 

540.2 

26.28 

0.00 

3.88 

93 

0.0002 

bd 

3 

117 

537.8 

25.93 

0.01 

4.27 

123 

0.0002 

bd 

3 

147 

546.0 

25.27 

0.01 

4.87 

153 

0.0002 

bd 

3 

178 

540.8 

24.74 

0.00 

5.24 

183 

0.0002 

bd 

3 

203 

530.9 

23.70 

0.00 

5.59 

208 

0.0003 

bd 

3 

228 

532.6 

22.87 

0.00 

6.02 

233 

0.0003 

bd 

3 

253 

533.8 

22.50 

0.00 

6.13 

258 

0.0003 

bd 

3 

278 

548.7 

22.45 

0.00 

283 

0.0003 

bd 

3 

303 

543.9 

21.58 

0.00 

6.96 

308 

0.0003 

bd 

3 

328 

558.1 

23.97 

0.00 

6.69 

333 

0.0004 

bd 

4 

10 

564.2 

29.27 

0.00 

3.07 

15 

0.0002 

bd 

443 

4 

50 

552.7 

27.98 

0.00 

2.85 

55 

0.0006 

bd 

4 

80 

544.8 

27.35 

0.00 

3.12 

85 

0.0009 

bd 

4 

105 

542.3 

27.32 

0.02 

3.32 

110 

0.0012 

bd 

4 

130 

539.7 

27.15 

0.07 

3.43 

135 

0.0016 

bd 

4 

155 

537.8 

26.20 

0.17 

3.68 

160 

0.0021 

bd 

4 

180 

560.7 

26.18 

0.43 

4.24 

185 

0.0028 

bd 

4 

205 

542.6 

23.85 

1.27 

5.29 

210 

0.0043 

bd 

4 

230 

548.8 

22.09 

2.74 

6.53 

235 

0.0054 

22.3 

241 

4 

255 

549.2 

19.44 

1.14 

7.94 

260 

0.0060 

28.8 

206 

4 

280 

549.6 

17.05 

3.01 

9.29 

285 

0.0088 

27.8 

315 

4 

305 

543.7 

14.18 

3.97 

10.83 

310 

0.0116 

55.2 

211 

4 

330 

544.2 

11.06 

5.39 

12.76 

335 

0.0183 

60.2 

303 

4 

355 

558.2 

9.21 

4.90 

13.64 

360 

0.0208 

71.9 

289 

4 

380 

542.6 

4.54 

5.33 

16.83 

385 

0.0362 

110.8 

326 

4 

405 

538.0 

1.59 

3.49 

19.68 

410 

0.1100 

122.5 

899 
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a 


Chloride 

(mM) 


Sediment 

Depth 

(cm) 


Core 

Pore 

Water 

Sed 

Depth 

(cm) 

Chloride 

(mM) 

Sulfate 

(mM) 

Sulfide 

(mM) 

DIC 

(mM) 

Gas 

Sample 

Sediment 

Depth 

(cm) 

ch4 

(mM) 

c2h6 

(nM) 

c,/c2 

(vol) 

SMI 

(cm) 

31 

265 

560.3 

22.85 

2.14 

6.59 

270 

0.00627 

bd 

31 

290 

544.3 

21.57 

2.59 

8.14 

295 

0.0 

bd 

32 

5 

544.7 

28.16 

0.01 

2.63 

10 

0.0 

bd 

654 

32 

35 

536.0 

28.07 

0.03 

2.67 

40 

0.0 

bd 

32 

60 

541.1 

28.57 

0.02 

2.61 

65 

0.0 

bd 

32 

85 

541.0 

28.62 

0.00 

2.70 

90 

0.0 

bd 

32 

110 

542.5 

28.60 

0.03 

2.75 

115 

0.0 

bd 

32 

135 

541.9 

28.62 

0.04 

2.77 

140 

0.0 

bd 

32 

160 

539.7 

28.57 

0.02 

2.87 

165 

0.0 

bd 

32 

185 

540.6 

28.50 

0.05 

2.83 

190 

0.0 

bd 

32 

210 

540.6 

27.77 

0.06 

3.22 

215 

0.0 

bd 

32 

235 

545.2 

27.57 

0.10 

3.63 

240 

0.0 

bd 

32 

260 

545.6 

26.36 

0.33 

4.76 

265 

0.0 

bd 

32 

285 

538.8 

25.51 

1.49 

5.41 

290.00 

0.00 

bd 

32 

310 

545.7 

22.94 

3.09 

6.99 

315.00 

0.01 

bd 

32 

335 

540.7 

19.90 

2.77 

8.25 

340 

0.0 

bd 

32 

360 

549.4 

19.39 

3.41 

9.36 

365 

0.0 

bd 

32 

385 

537.0 

17.93 

4.02 

9.85 

390 

0.0 

bd 
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Appendix  II.  Coring  Notes 


Core  sample  distributions  were  outlined  for  preliminary  discussions  on  the  sample  processing. 


10  cm 


Pore  water  parameters 

Volume 

Cl,  so4'2 

2 

h2s 

1 

DIC 

2 

813C-DIC/DOC/613DOC/ 

Flourescence 

3 

ICP 

3 

TA 

1 

Organic  Acid 

2 

Nutrients 

7 

Br,  I 

4 

Archived  Samples 

TOTAL 

25 

top  section 


Sediment  parameter 

Mass  gm 

Carbonate 

%  Organic  C,  TCN 

10 

Community  Diversity, 

sampled  prior  to  press 

10 

[CnHn],  513C,  sampled  prior  to 

press 

3 

Porosity,  sampled  prior  to 

press 

1 

Hopanes 

10 

10  cm 


Pore  water  parameters 

Volume 

no  pore  water  providing 
volume  above  is  sufficient 

bottom  section 


Sediment  parameter 

Amount 

Microbiology 

Biomarkers 

5  cm  section 

This  figure  represents  one  section  of  a  core  for  the  sample  distribution 
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The  preliminary  core  log  sheet  changes  through  the  cruise  according  to  core  lengths. 


Cruise 

Core  ID 

Lat 

Long 

Date 

Time 


Comments 


Observer 


Top  Top  Top 
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Appendix  III  Cruise  Schedule 

Times  in  hours,  total  maximum  time  available  12.5  days  at  sea  =  300  hrs. 

1 .  Transit  Wellington/Study  Area  A:  Time  Total 
To  Waypoint  cw-1 

-130  nm,  +  1  h  in  port,  10  knots  14  14.0 

2,  “Warm-up”  core  and  CTD  deployment 
Study  Area  A,  Objectives  2,  5,  7,  Priority  1 
Approx.  600-800  m  water  depth 


3.5  kHz,  12  kHz,  and  multibeam  Waypoints  cw-1  -  cw-2 

8.1  nm  at  10  knots  0.8 

Back  to  cw-3  0.4 

CTD  at  Waypoint  cw-3  1 .5 

Change  to  coring  0.5 

Gravity/piston  core  Waypoint  cw-3  or  based  on  3.5  kHz  data  1 .5 
Seismic  survey  starts  close  to  cw-3 

(sub-total)  4.7  18.7 


3.  Seismic  survey  A1 

Study  Area  A,  Objective  1,  Priority  1 

45/105  cu-in  GI  gun,  c.  25  m  shot  spacing,  0.5  msec  sampling 

Source/streamer  depth  5  m 

Also  run  3.5  kHz,  12  kHz,  and  multibeam 


Deploy  gear  4 

Survey  Waypoints  sal-1  -  sal-8 

89  nm  at  4.5  knots  19.8 

6  turns  at  10  additional  minutes  1 

Retrieve  gear  1 

(sub-total  25.8)  44.5 


4.  Coring,  heatflow,  and  CTD  transect  A1 
Study  Area  A,  Objectives  2-5  Priority  1 
Transect(s)  between  Waypoints  cal-1  and  cal -2 
-2000-2500  m  water  depth 

Note:  coring  includes  “failed”  cores;  heatflow  and  CTD  numbers 
are  essentially  placeholders;  coring  will  probably  be 
broken  up  into  at  least  2-3  segments 


Coring  11  Stations  (3  hrs.  each),  10  transits  (0.5  hrs.)  38 

Heatflow  6  Stations  (3  hrs.  each),  5  transits  (0.5  hrs.)  20.5 

CTD  8 

(sub-total)  66.5)  111.0 

5.  Transit  Study  Area  A  /  Temperature  sensor 
To  Waypoint  t_probe 

80  nm  at  10  knots  8  1 19.0 
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6.  Recover  temperature  sensor 
Study  Area  B,  Objective  6,  Priority  1 
Waypoint  t_probe,  -600  m  water  depth 


2 


121.0 


7.  Transit  Paleo-core  1  stations 

Note:  Location  can  change  significantly 
To  Waypoint  cpl-1 

7.9  nm  at  10  knots  0.8  121.8 

8.  Paleo-core  1 

Note:  Location  can  change  significantly 
Study  Area  B,  Objective  7,  Priority  1 
-800  m  water  depth 


3.5  kHz,  12  kHz,  and  multibeam  Waypoints  cpl-1  -  cpl-2 

13.5  nm  at  10  knots 

1.4 

Coring  at  Waypoint  cpl-3  and/or  additional  based  on  3.5  kHz 
Coring  3  Stations  (2  hrs.  each)  2  transit  (0.5  hrs) 

7 

(sub-total 

8.4) 

130.2 

9.  Transit  to  BOL  of  Seismic  Survev  B2 

To  Waypoint  sb2-l 

17.6  nm  at  10  knots 

1.8 

132.0 

10.  Seismic  survev  B2 

Note:  only  if  sufficient  time  to  complete  remaining 

Priority  1  experiments  (dr edging^ coring  below) 

Study  Area  B,  Objective  9,  Priority  2 

GNS’  Mini-GI  gun  (chamber  configuration  tbd) 

-25  m  shot  spacing,  0.5  msec  sampling 

Source/streamer  depth  ideally  2.5  m,  but  may  need  to  be 

deeper  unless  seas  are  unusually  calm. . . 

Deploy  gear 

Survey  Waypoints  sb2-l  -  sb2-19 

4 

84.0  nm  at  4.5  knots 

18.7 

18  turns  at  10  additional  minutes 

3 

Retrieve  gear 

Dredging  close  to  EOL 

1 

(sub-total 

26.7) 

158.7 

1 1 .  Dredging  and  coring.  Priority  1 

Note:  Conducted  before  seismic  survey  above  if  time  is  tight 
Study  Area  B,  Objective  8,  Priority  1 

Dredge  at  Waypoint  dbl  -600  m  water  depth 

3 

Transit  to  core  location  cbl,  7.8  nm  at  10  knots 

0.8 

Coring  cbl  -2700  m  water  depth 

3 

(sub-total 

6.8) 

165.5 

38 


(End  of  Priority  1 — note:  w/o  Seismic  Survey  B2 

Transit  back  to  Wellington 

->  140.2  hours  contingency  ->  47%  downtime) 
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138.8 

159.8 


12.  Dredging  and  Coring  Priority  2 
Study  Area  B,  Objective  8,  Priority  2 


Transit  from  Waypoint  cb  1  to  db2,  11.32  nm  at  1 0  knots  1 . 1 

Dredge  at  Waypoint  db2  -600  m  water  depth  3 

Transit  to  Waypoint  cb2,  2.3  nm  at  10  knots  0.2 

Coring  cb2 -1400  m  water  depth  2 

(sub-total  6.3)  171.8 

13.  Transit  Study  Area  B 
To  Waypoint  cal-2 

60  nm  at  10  knots  6  177.8 


14.  Coring,  heatflow.  and  CTD  transect  A2 
Study  Area  A,  Objectives  2-5  Priority  2 

Locations  tbd,  likely  transect(s)  between  Waypoints  cwl  and  cal-2 
-500-2500  m  water  depth 

Note:  1.  coring  includes  “failed”  cores;  heatflow  and  CTD  numbers 
are  essentially  placeholders;  coring  will  probably  be 
broken  up  into  at  least  2-3  segments 

2.  Can  be  shortened  to  allow  further  Priority  2  work 

3.  Paleo-core  2  may  be  done  this  transect 


Coring  1 1  Stations  (3  hrs.  each),  10  transits  (0.5  hrs.)  38 

Heatflow  6  Stations  (3  hrs.  each),  5  transits  (0.5  hrs.)  20.5 

CTD  8 

(sub-total)  66.5)  244.3 


15.  Paleo-core  2 
Note:  Location  tbd 

May  be  cored  before  Pore-water  coring/heatflow  transect 


Study  Area  A?,  Objective  7,  Priority  1 
-800  m  water  depth? 

Transit  to  site  (placeholder)  2 

3.5  kHz,  12  kHz,  and  multibeam  1.5 

Coring  3  Stations  (2  hrs.  each)  2  transit  (0.5  hrs)  7 

(sub-total  10.5)  254.6 

(End  of  Priority  2  —  use  of  any  remaining  time  tbd, 
likely  in  Study  Area  A) 

1.  Transit  Study  Area  A/Wellington 

(see  1)  14  268.6 


->31.4  hrs  for  300  contingency  (10.5%) 


39 


Appendix  IV.  On  Board  Instrumentation. 

NIW  A/Vessel: 

•  SIMRAD  EM  300  multibeam  echosounder 

•  SIMRAD  EA  500  echosounder  (12  kHz,  38,  120  kHz),  logging  software  by  GNS,  see  below) 

•  SIMRAD  EA  60  echosouner  (on  bridge) 

•  Knudsen  3.5  kHz  subbottom  profiler 

•  Standard  piston/drop  corer 

•  Multi-channel  seismic  system,  GI  gun  (total  210  cu-in,  preferred  configuration  45/105  cu-in  in  GI 
mode),  Hydroscience  oil-filled  streamer  (500m,  6.25  m  hydrophone  spacing,  12.5  m  group 
spacing),  Geometries  recoder  (preferred  0.5  msec  sampling  rate) 

•  CTD  system 

•  Laboratory  containers  (2)  20-ft 
GNS 

•  GNS  dredge  (stored  at  NIWA) 

•  GNS  mini-GI  gun  with  spares  (chamber  configuration  tbd) 

•  Gas  chromatography  for  CH4  analyses  (including  H-generator  and  laptop  for  data  recording) 

•  Wet  geochemical  equipment  for  nutrients  (including  photometer,  titrator  unite,  lab  ware) 

•  Linux  PC's  for  data  logging  and  processing  including  Claritas  seismic  software  (GNS  private 
network) 

•  Plotter  (GNS  network) 

•  Windows  PC  for  logging  EK500  data  with  Echoview  (ship’s  network) 

NRL/U  Hawaii/U  of  Rochester 

•  Heatflow  probe 

•  GC  for  CH4  to  C3H8 

•  Pore  water  press 

•  Ion  chromatography  for  pore  water  analyses  (including  lab  ware) 

•  Equipment  for  microbiological  sampling 

•  Fluorometer  for  sulfide 

•  Coulometer  for  dissolved  inorganic  carbon 

•  Deck  and  lab  tools 

Otago 

•  Linux  laptops  for  data  logging  and  processing  (to  hook  into  GNS  network) 
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Appendix  V  On  Board  Research  Team  Breakout 


NRL  coring/heatflow  team: 

PI:  Coffin 

Scientists:  Wood,  Hamdan,  Pohlman,  Smith,  Dale,  MacDonald,  Downer, 
Bryant,  Kurasaki,  Yoza 

GNS  coring: 

PI:  Greinert 

Scientists:  Schinteie,  Tomaru,  (+Faure,  Kurazaki,  Yoza,  Kukowski, 
if  available) 

Seismic  and  3.5  kHz: 

(in  addition  to  operation  conducted  by  NIWA  personnel) 

PI:  Gorman 

Scientists  Crutchley,  Senger,  Kilner,  Fohrmann,  (+Wood,  Pecher,  Henrys, 
if  available) 

(beyond  here,  very  flexible) 

GNS  water-column  chemistry  (CTD): 

PI:  Faure 

Scientists:  (Greinert,  if  available) 

Paleo-coring: 

PI:  Orpin 

Scientists:  tbd  (+Schinteie,  Kurazaki,  Yoza,  Tomaru,  Kukowski,  if  available) 
12  kHz: 

(in  addition  to  operation  conducted  by  NIWA  personnel) 

PI:  Greinert 

Scientists:  tbd,  if  needed 
Multibeam: 

(in  addition  to  operation  conducted  by  NIWA  personnel) 

PI:  Kukowski 

Scientists:  tbd,  if  needed  (+Orpin,  if  available) 
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Appendix  VI.  NIWA  Data  Log  Sheet 
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